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ABSTRACT 

We study the observational constraints of CMB temperature and polarization anisotropics 
on models of dark energy, with special focus on models with variation in properties of dark 
energy with time. We demonstrate that the key constraint from CMB observations arises from 
the location of acoustic peaks. An additional constraint arises from the limits on Qnr from the 
relative amplitudes of acoustic peaks. Further, we show that the distance to the last scattering 
surface is not how the CMB observations constrain the combination of parameters for models 
of dark energy. We also use constraints from Supernova observations and show that unlike 
the Gold and Silver samples, the SNLS sample prefers a region of parameter space that has a 
significant overlap with the region preferred by t he CMB observations. This is a verific ation 
of a conjecture made by us in an earlier work (Jassal, Bagla, & Padmanabhan 2005a). We 
discuss combined constraints from WMAP5 and SNLS observations. We find that models 
with w ~ — 1 are preferred for models with a constant equation of state parameters. In case 
of models with a time varying dark energy, we show that constraints on evolution of dark 
energy density are almost independent of the type of variation assumed for the equation of 
state parameter. This makes it easy to get approximate constraints from CMB observations on 
arbitrary models of dark energy. Constraints on models with a time varying dark energy are 
predominantly due to CMB observations, with Supernova constraints playing only a marginal 
role. 
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confirmation usi ng o b servations of high r e dshift supernovae 



i Garnavich et al. 



I Sahni & Starobinskvl 120001; iPadmanabhan l 120031; iPeebles & Rata 
20031; lEllij 120031; IPadmanabhanl 12005 j ; Iperivolaropoulosl |200: 
ICopeland. Sami & Tsuiikawal 1 20061) and it is known to be 
consistent with observations. In order to avoid theoretical 
problems related to c o smolo gical constant ( [Weinberg! 1 19891 ; 
ICarroll. Press, & Turner! 1 19921) . many other scenarios have 
been investigated: these i n clude quintessence ( Steinhardt 



— — i i ir- 1 i 1 — ii 1 oeen investigates tnese nicitiue quintessence < otemiiaiu 

1998; Perlmutter et al. 1999; Tonrv et al. 2003; dg k & Urena . Lopez & M atos 



iBarris et alj|2004l ; iRiess et alj|2004l ; lAstier et alj|2005h has made 



2003 



2000 



this result more acceptable to the community. Using these obser- 
vations along with observations of cosmic microwave background 
radiation (CMB) jMelchiorri et al.l |2000| ; ISpergel et all 120031; 
Komatsu et all [2009) and large scale structure iPercival et al.l 
we can construct a "concordance" model for cos mology 



Gonzalez-Diaz 2002; de Ritis & Marino 2001; Sen & Seshadri 
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and study var i ations around it (e . g., see Dunklev et al. (2008); 
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for an overview of our current understanding, see iPadmanabhai 



34); 
han 



iKomatsu et all d2009h : iBridle et al] O003) ; Tegmark et al. (2004); 
for an overv 
d2005rffl ) 

Observations indicate that the dominant component of energy 
density — called dark energy — should have an equation of state 
parameter w = P/p < —1/3 for the universe to undergo acceler- 
ated expansion. Indeed, present day observations require w ~ — 1. 
The cosmological constant is the simplest explanation for accel- 
erated expansion dWeinberdll98^ : ICarroll. Press. & Turned fl992l ; 



JArmendariz-Picon. Mukhanov, & Steinhardi 



12003: iMakmarti et all 120031; IChimento & Feinstein 
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Sen & Scherreil 2005 ), ph antom fields 1 Caldwell 
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Onemli & Woodard 2002; Carroll. Hoffman & Trodden 



Singh. Sami & Dadhicb 120031 ; iFramptonl 



2003; 



Jassal, Bagla and Padmanabhan 



Gonzalez-Diadl2003hlDabrowski. Stachowiak & Szvdlowskil2003l; 
Elizalde. Noiiri. & Odintsov 20041; iNoiiri. Odintsov. & Tsuiikawe 



2005 



2007 



2001 
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Bronnikov & Starobinsky 2007), branes dUzawa & Soda 
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2005J12002I; 



2004 



Sahni & Shtanovl I2OO3I) 




Onemli & Wood arcj 120041; 



Andrianov. Cannata & Kamen shchik 
Lazkoz. Neseseris & Perivolaropoulosl 20051; 



20051; 



ZUUa: ILaz Koz. INesesens & Perivolaropoulosl UUUa; 

Cognola. Elizalde. Noiiri. Odinstov & Zerbmj 20061; iRen & Mend 



Polarski & Ranquet 2005; 
Das, Banerjee & Dad hichl 120051; 
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Sola & Stefancid 120051; 
Apostolopoulos & Tetradisl 



Arianto et al.l 120071 ; ISahni. Shafieloo & Starobinskvl 120081) 
In these models one can have -1 and in general w varies with 
redshift. For references to papers that discus s specific models, the 
reade r may consult one of the many reviews dSahni & Starobinskyl 
200d; IPadmanabhanl [20031; IPeebles & Ratral 120031 ; lEllij 120031 ; 
Padmanabhan! l2005al ; lAlcanizll2006l) . Even though these models 



have been proposed to overcome the fine tuning problem for 
cosmological constant, most models require similar fine tuning 
of parameter(s) to be consistent with observations. Nevertheless, 
they raise the possibility of w(z) evolving with time (or it being 
different from —1), which can be tested by observations. 

Given that w for dark energy should be smaller than — 1 /3 for 
the Universe to undergo accelerated expansion, the energy density 
of this component changes at a much slower rate than that of matter 
and radiation. (Indeed, w = — 1 for cosmological constant and in 
this case the energy density is a constant.) Unless w is a rapidly 
varying function of redshift and becomes w ~ at (z ~ 1), the 
energy density of dark energy should be negligible at high redshifts 
(z S> 1) compared to that of non-relativistic matter. If dark energy 
evolves in a manner such that its energy density is comparable to, or 
greater than the matter density in the universe at high redshifts then 
the basic structure of the cosmological model needs to be modified. 
We do not consider such models, instead we confine our attention 
to constraints on dark energy in realistic models and choose obser- 
vations which are sensitive to evolution of w(z) at redshifts z < 1. 

Supern ova observations pe rmit a large variation in the equa- 
tion of state dAlam eta l. 2004a. b). It has recently been argued that 
a cosmological constant with a small curvature can be interpreted 
as a d ynamical dark en ergy model dClarkson. Cortes. & Bassettl 
120071 ; I Virev et all 120051) . We have shown that a combination 
of supernova observations with CMB observations and abun- 
dance of rich clusters of g alaxies provides tight constraints on 
variation of dark energy dJassal. Bagla & Padmanabhanl 120051 ; 
Ijassal. Bagla. & Padmanabhan! l2005al) . Of these CMB data pro- 
vides the most stringent bounds on the allowed variation in evo- 
lution of dark energy density. In this paper the main motivation is 
to study how these models fare in the lig ht of current CMB data 
jPunklev et al.ll2008l ; lKomatsu et alj2009t) . 

A variety of observations ca n be used to constrain models 
of dark energy, e.g. see §11 B of Ijassal. Bagla. & Padma nabhanl 
(2005a) for an overview. Observations of high redshift su- 
pernovae provided the firs t direct ev i dence for accelerated 
expan sion of the universe dRiess et ail 1 19981 ; IPerlmutter et al.l 
1 19991) . This, coupled with the ease with which the high 
redshift supernova data can be compared with cosmological 
models has made it the favourite benchmark for compari- 
son with models of dark energy. It is often considered suf- 
ficient to compare a model with the supernova data even 
though observers a nd theorists have pointed out potential prob- 
lems with the data dPerlmutter and Sch midt 20031; Pain & Ralstonl 



120051) as well as some peculia r implications of the data 
(Padmanabhan & ChoudhurJ 120031 ; Ichoudhurv & Padmanabhanl 
120051; IShashikantll2005l) . Further, it has been shown that other ob- 
servations like temperature anisotropies in the CMB fix the distance 
to the last scattering surface and are a reliable probe of dark ene rgy 
dEisenstein & Whitel2004l ; |jassal, Bagla & Padmanabhanl2005h . In 
this work, we try to understand the origin of the constraint on mod- 
els of dark energy from CMB anisotropy observations. We show 
that the angular scale of acoustic peaks provides the leading con- 
straint on the combination of parameters. An additional constraint 
comes in through the degeneracy with £Inr, where a constraint on 
Q nr from the relative heights of acoustic peaks translates into an 
indirect constraint on the equation of state parameter w. 

In an earlier work, we compared the constraints on models of 
dark energy from supernova and CMB observations and pointed 
out that models preferred by these observations lie in distinct parts 
of the parameter space an d there is no overlap of regions allowe d 
at 68% confidence level dJassal, Bagla. & P admanab hanl l2005ah . 
Even though different observational sets are sensitive to different 
combinations of cosmological parameters, we do not expect models 
favoured by one observation to be ruled out by another when such 
a divergence is not expected. This divergence may point to some 
shortcomings in the model, or to systematic errors in observations, 
or even to an incorrect choice of priors. We suggested that this may 
indicate unresolved systematic errors in one of the observations, 
with supernova observations being more likely to suffer from this 
problem due to the heterogeneous nature of the da ta sets available a t 
the time. In Supernova Legacy Survey (SNLS) dAstier et alj|2005l) 
survey, a concerted effort has been made to reduce systematic errors 
by using only high quality observations. The systematic uncertain- 
ties are reduced by using a single instrument to observe the fields. 
Using a rolling search technique ensures tha t sources are not lo st 
and data is of superior quality (for details see lAstier et al.l d2005)). 
If our claim about Gold+Silver data set were to be true, SNLS data 
should not be at variance with the WMAP dattQ In this work we 
study constraints on dark energy models from high redshift super- 
nova observations from the SNLS survey and also observations of 
the temperature and polarisation anisotropies in the CMB using the 
WMAP5 data. 

This paper is a revised version of an earlier manuscript which 
became out of date after WMAP-3 and then WMAP-5 data were 
released. In the interim period, the issue of systematics in inhomo- 
geneous data sets has been accepted, therefore we do not emphasize 
that aspect much in this version. The focus of the current paper is 
twofold: to study constraints on dark energy models in light of the 
SNLS and WMAP-5 data, and, to understand the combination of 
cosmological parameters that is constrained by the CMB observa- 
tions. 



2 DARK ENERGY 

2.1 Cosmological equations 

If we assume that each of the constituents of the homogeneous and 
isotropic universe can be considered to be an ideal fluid, and that 
the space is flat, the Friedman equations can be written as: 



(:)' 



8ttG 



(1) 



1 This has been shown by several authors, including ourselves in a much 
earlier version of this manuscript. 
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| = - 4 -fiP + 3P) (2) 

where P is the pressure and p — pnr + p~ t + p DE with the re- 
spective terms denoting energy densities for nonrelativistic matter, 
for radiation/relativistic matter and for dark energy. Pressure is zero 
for the non-relativistic component, whereas radiation and relativis- 
tic matter have P 7 = p 7 /3. If the cosmological constant is the 
source of acceleration then p DE — constant and P DE = — p DE . 
Analysis of nonflat cosmologies reveals that allowed range of cur- 
vature of the universe Qk is —0.012 — 0.009 at 95% c onfidence 
level and a flat universe is a good fit to the current data dXia et all 
l2008h . 

An obvious generalisation is to consider models with a con- 
stant equation of state parameter w = P/p = constant. One can, in 
fact, further generalise to models with a varying equation of state 
parameter w(z). Since a function is equivalent to an infinite set of 
numbers (defined e.g. by a Taylor-Laurent series coefficients), it is 
clearly not possible to constrain the form of an arbitrary function 
w(z) using a finite number of observations. One possible way of 
circumventing this issue is to parameterise the function w(z) by a 
finite number of parameters and try to constrain these parameters 
with the available observational data. There have been many 
attempts to describe varying dark energy with different parame- 
terisations ( Wang & Tegmark 2004; Bassett, Corasaniti & Kunz 
20041; IJassal, Bagla & Padmanabhanl |2005l: IL"eel [20051: 10 120041: 
Hanne stad & Mortselll |2004|; Ijassal. Bagla. & Padmanabhanl 
2005a) where the functional form of w{z) is fixed and the variation 
is described with a small number of parameters. Observational 
constraints depend on the specific parameterisation chosen, but it 
should be possible to glean some parameterisation independent 
results from the analysis. 

To model varying dark energy we use two parameterisations 

w(z) = w + w'(z = 0) T Z ; p = l, 2 (3) 
(1 + z) p 

These are chosen so that, among other things, the high redshift 
behaviour is completely diffe r ent in these two parameterisations 
Jassal, Bag la & Padmanabhanl {2005). If p = 1, the asymptotic 
value w(oo) — wo + w'(z — 0) and for p = 2, w(oo) = Wo, 
For both p — 1,2, the present value w(0) = Wo. Clearly, we must 
have w(z 3> 1) < —1/3 for the standard cosmological models 
with a hot big bang to be valid. This restriction is imposed over and 
above the priors used in our study. 

The allowed range of parameters wq and w' = w'(z = 0) 
is likely to be different for different p. However, the allowed vari- 
ation at low redshifts in po e should be similar in both models as 
observations actually probe the variation of dark energy density. 

2.2 Observational Constraints 

2.2.1 Supernova Data 

In this work, we concentrate on SN and WMAP observations. 
SN data provides geometric constraints for dark energy evolution. 
These constraints are obtained by comparing the predicted lumi- 
nosity distance to the SN with the observed one. The theoretical 
model and observations are compared for luminosity measured in 
magnitudes: 

m B (z) = M + 5log w (D L ) (4) 

where M = M — 5logio(Ho) and Dl = HocIl, M being the 
absolute magnitude of the object and <1l is the luminosity distance 
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w 

Figure 1. This figure shows contours of angular diameter distance. 
The range plotted is the range allowed by WMAP5 data. From the 
line on to the one at bottom the values correspond to cIa = 
12000, 13000, 14273, 20000 and 20100 Mpc. 
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Figure 2. This figure shows contours of constant angular size of the 
Hubble radius. The contours from left to right correspond to 9~ 1 = 
0.0165, 0.017, 0.18, 0.020. 



d L = (1 + z)a(t )r(z); r(z) = cj -^y (5) 

where z is the redshift. This depends on evolution of dark en- 
erg y through H(z). For our analysis we use the SNLS data 
set jAstieretal.ll2005h and for reference we also used the com- 
bined gold and silver SN data set in iRiess et aL I d2004h (see also 
dNesseris & Perivolaropoulosll2005h). Thi s data i s a collection of 
supern ova observations from iTonrv et alj d2003l) : iGarnavich et al.l 
dl998h and many other s ources with 16 su pernovae discovered with 
Hubble space telescope dRiess et al.l2004l) . The parameter space for 
comparison of models with SN observations is small and we do a 
dense sampling of the parameter space. 

2.2.2 CMB Data 

CMB anisotropies constrain dark energy in two ways, through the 
distance to the last scatte ring surface and throug h the Integrated 
Sachs Wolfs (ISW) effect dPeiris and Spergell2000h . Given that the 
physics of recombination and evolution of perturbations does not 
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Figure 3. Marginalized likelihood contours for different parameters for ACDM model. The regions enclosed by the contours are 68% and 95% confidence 
limits. The results are consistent with WMAP5 results for ACDM model. The left plot shows the allowed range in n — Q^r R h 2 plane. The next figure shows 
the correlation between parameters n and £lgh 2 . The figure on the right shows corresponding contours inn — t plane. 



change if w(z) remains within some safe limi ts, any change in the 
locati on of peaks will be due to dark energy ( lEisenstein & White! 
|2004) . For models with a vari able w(z), the constraint is essentiall y 
on an effective value u> e // (Jassal, Bagla, & Padmanabhan 2005a). 
This constraint can arise either through the angular diameter dis- 
tance or the angular size of the acoustic horizon seen reflected in 
the scale corresponding to the acoustic peaks in the CMB angular 
power spectrum. 

Figure 1 shows contours of equal angular diameter distance 
to the last scattering surface in the Q. nr — w plane. We assumed a 
fixed value of Ho and Qb for this plot. If the distance to the last 
scattering surface is the key constraint on models of dark energy 
then the likelihood contours should run along the contours of con- 
stant distance in this plane. 

We may use the angular size of the Hubble radius at the time 
of decoupling as an approximate proxy for the angular size of the 
acoustic horizon for the purpose of this discussion. The approxi- 
mate angular size 9 of the Hubble radius at the time of decoupling 
can be written as: 

H /H(z) 



fdy/(H(y)/H ) 



(n NR {i + z f) 



3N-1/2 



fdy/y/{l NR (l + zf + g°z( Z )/Q°i 



(n NR (i + z ) 3 y 1/2 



j d y /Jn NR (i + zf + n de (i + z f(i+™ef f ) 



(6) 



Clearly, the value of the integral will be different if we change wo, 
w'(z — 0) and there will also be some dependence on the parame- 
terized form. If the location of peaks in the angular power spectrum 
of the CMB provide the main constraint, this can only constrain 
W e ff and not all of wq, w'(z — 0) and p. Therefore if the present 
value wq < w e f f then it is essential that w'(z — 0) > 0, and sim- 
ilarly if wo > w e f f then w'(z = 0) < is needed to ensure that 
the integrals match. Specifically, the combination of wq, w'(z = 0) 
and p should give us to e // within the allowed range. 

Figure 2 shows contours of constant 9 in the Q,nr — w plane. 



If the angular scale of acoustic peaks is the key constraint arising 
from CMB observations for models of dark energy then we should 
see two features in the likelihood contours: 

• Likelihood contours for models with a constant equation of 
state parameter w in the Q, mr — w plane should run along contours 
of constant 9, as shown in Figure 2. 

• Likelihood contours for different models of varying dark en- 
ergy should coincide in the Q.nr — w e ff plane. These should also 
coincide with the contours for models with a constant w in the 
fijvi? — w plane. 

The origin of the CMB constraints on dark energy therefore is not 
in the raw distance to the surface of last scattering but in the com- 
bination of parameters that determines the location of peaks in the 
angular power spectrum. It is important to note that the distance to 
the surface of last scattering is a derived quantity. The CMB obser- 
vations constrain only one number, the effective equation of state. 
There is no ambiguity e xcept for models with early dark energy 
dLinder & Robbersll2008l) . In these models, the growth of perturba- 
tions is slo wer than models in which dark energy comes into play a t 
late times teenabed & Bernardeaj|200ll : lDoran & Robbers| [2006). 

In our analysis, we use the ang ular power spectrum 
of the CMB temperature anisotropies dHu & Dodelsonl 1 20021 : 
I White & Cohnll2002l:[Subramanianll2004) a s observed by WMAP 
jDunklev et"ai] 120081 : licomatsu et alj |2009» and these are com- 
pared to theoretical predictio ns using the likelihood program pro- 
vided by the WMAP team jDunklev etai] 120081 : iKomatsu etail 
2009). We vary the amplitude of the spectrum till we get the 
best fit with WMAP obs ervations. The CMBFASlQ package 
(Seljak& Zaldarriaga 1996) is used for computing the theoretical 
angular power spectrum for a given set of cosmological parame- 
ters. We have combined the likelihood program with the CMB- 
FAST code and this required a few minor changes in the CMB- 
FAST driver routine. We also made changes in the driver program 
to implement Monte C arlo Markov Chain for sampling the p aram- 
eter space. Please see Jjassal, Bagla, & Padmanabhan] l2005ij) for 
details of the MCMC implementation. 

Although we can use other observations like abundance of 
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Figure 4. This figure shows contours of angular diameter distance 
(blue/dashed) and of constant angular size of the Hubble radius (red/solid) 
overlaid with likelihood regions allowed by WMAP5 data in the Qnr ~ w 
plane. 



Table 1. This table lists the priors used in the present work. Apart from 
the range of parameters listed in the table, we assumed that the universe is 
flat. We assumed that the primordial power spectrum had a constant index. 
Further, we ignored the effect of tensor perturbations. The range of values 
for wq and w'(z = 0) is as given below, but with the constraint that w(z = 
1000) < — 1/3. Any combination of wq and w'(z = 0) that did not satisfy 
this constraint was not considered. Values given in parenthesis were used for 
analysing constraints from high redshift supernovae in §3.1. 



Parameter 


Lower limit 


Upper limit 


n B 


0.03 


0.06 




0.1 


0.5(0.7) 


h 


0.6 


0.8 


T 


0.0 


0.4 


n 


0.86 


1.10 


w 


-2.0(-4.0) 


-0.4 


w'{z = 0) 


-5.0 


5.0 



rich clusters, baryonic features in the power spectrum, etc. but we 
find that the two observations used here are sufficient for this study 
(Jassal 2009). 



3 RESULTS 

In this section we will describe the results of our study. We studied 
models in three classes: 

• Models with a constant equation of state parame- 
ter w . We studied mode l s with perturbations in da r k en- 
ergy dBean & Dorel |2004l: ICaldwell. Dave & Steinhardtl 1 19981: 
Welle r & Lewis! 12003: iHannestaol 120051: iFabris. Shapiro. & Sola 



20071 : Fabris & Goncalvedl200d: lUnnikrishnan. Jassal. & Seshadril 
20081: IJassall 120091: iBartolo. Corasaniti. Liddle & MalquartfeOolF 
Mota & van de Brucki2004l:lGordon & Hul2004l : |Gordon & Wandsl 
20051: iNunes :& Mot A l2006l : ISergiienko et alj 12008 
Maininill2009[ 2008) as well as without. 
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• Models with a varying equation of state parameter w, with 
variation given by Eqn.3 (p=l). Perturbations in dark energy were 
not taken into account in this case. 

• Models with a varying equation of state parameter w, with 
variation given by Eqn.3 (p=2). Perturbations in dark energy were 
not taken into account in this case too. 

We analyse the allowed range of cosmological parameters for 
these cosmologies and consider the probability with which the 
ACDM model is allowed within these three classes of mod- 
els. In light of the significant disagreement between the al- 
lowed range of parameters from the high redshift supernova 
data from the Gold+Sil ver set and the CMB anisotropics 
from WMAP observations Jjassal. Ba gla. & Padmanabhan 2005a; 
iNesseris & Perivolaropoulosl 20051) . we also check the degree of 
overlap between the parameter space allowed by the supernova and 
the CMB observations respectively. The newly released 'Union' 
dataset includes d ata from the Supernova Legacy Survey, the 
ESSENCE Survey dWood-Vasev et alj|2007l) and the extended dis- 
tant supernova datas et from HST along with the older datasets 
dKowalski etai72 008). The combined data favours variation in dark 
energy equation of state. 



3.1 Cosmological Constant 

We begin with a very brief review of constraints on parameters in 
the case where a cosmological constant is the source of acceler- 
ating expansion of the universe. We used priors given in Table 1, 
except that the equation of state parameter is fixed to w — —1. 
Constraints on cosmological parameters are listed in Table 2. We 
would l ike to note that our results match those o btained by other 
authors (Dunkley et al. 20081: iKomatsu et aL .2009), Figure 3 shows 
contours of likelihood for some pairs of parameters as an example. 
We have shown contours in the n — jv nh 2 , n — Qsh 2 , and, n — r 
plane, we see that there is a strong correlation between n — Qsh 2 . 
These likelihood contours have also been shown as a reference for 
equivalent plots for models with t»/ -1, and help in checking the 
effect of the additional dark energy parameters on allowed range of 
other parameters. 



3.2 Constant w 

We first evaluate the nature of the CMB constraint on models of 
dark energy. Figure 4 shows likelihood contours from WMAP5 ob- 
servations in the £Inr — w plane for constant w models. We find 
that the orientation of these contours is roughly along contours of 
constant 6. To illustrate this, we have overlaid contours from Fig- 
ure 1 and Figure 2 in Figure 4. On the other hand there is no sim- 
ilarity between the likelihood and contours of distance to the last 
scattering surface. Thus we may conclude that the dominant con- 
straint provided by the CMB observations arises from the location 
of peaks in the angular power spectrum. The reason for this is that 
the angular diameter distance to the last scattering surface is a de- 
rived quantity, whereas the location of peaks in the angular power 
spectrum of temperature anisotropies is a direct observable. 

We use priors given in Table 1 for models with a con- 
stant equation of state parameter, with the obvious constraint that 
w'(z = 0) = 0. For supernova observations, we used wider pri- 
ors for w and Qnr in order to illustrate the differences between 
the two data sets studied here. We begin with a brief summary of 
results for the Gold+Silver data set. The best fit model in this case 
is w = —1.99 and £Inr = 0.47. The allowed range for w at 95% 
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Figure 5. Marginalized likelihood contours in Qjv\r — w plane for different models. In the top figures the blue dashed lines correspond to 68%, 95% 
confidence level using the gold data. The red solid lines correspond to confidence levels in SNLS data. The black thick lines are marginalized confidence levels 
using WMAP5 data. The top left panel if for constant w models with a homogeneous dark energy and the one on the right is when we include dark energy 
perturbations. The bottom left plot is for models with p = 1 and the right plot is with p = 2. 



confidence limit for large priors is —3.73 < w < —1.25. The cor- 
responding range for the density parameter is 0.28 < Qnr < 
0.57. With SNLS data, the best fit model is w = -1.06 and 
Qnr = 0.29. The allowed range for w at 95% confidence limit 
is —2.36 < w < —0.74. The allowed range for the density param- 
eter is 0.11 < Qnr < 0.48. There is clearly a large shift in the 
allowed values of parameters. 

We illustrate this in Figure 5 (top-left panel) where we have 
plotted the regions allowed by the two data sets at 68% confidence 
levels in the w — Q nr plane. Dashed line shows the region allowed 
by the Gold+Silver data set and the solid line is for the SNLS data 
set. We can deduce the following from this figure: 

• The region allowed by these two data sets at 68% confidence 
level has some overlap, thus we may say that the two sets are con- 
sistent with each other. 

• The overlap is at Qnr > 0.36 and is thus at margins of what 
is allowed by other observations. 

• The ACDM model is ruled out at 68% confidence level by the 
Gold+Silver data set. 

• The best fit of each set is ruled out by the other data set at 
this confidence level. Indeed, the best fit of Gold+Silver data set is 
allowed by the SNLS data with a probability V = 12.65% while 
the best fit of the SNLS data set is allowed by the Gold+Silver data 
set with? 3 = 8.14%. 



This point is reiterated by the likelihoods of w and Qnr for these 
models in the same panel. 

The figure shows the large overlap between the likelihood 
curves corresponding to SNLS data and WMAP data where the 
Gold+Silver data clearly favours higher values of Qnr and more 
negative w. The phantom models are still allowed but the SNLS 
data as well as WMAP data show a preference for models close to 
a cosmological constant. 

For comparison, WMAP allows —1.25 < w < —0.7 and 
0.2 < Qnr < 0.38 if dark energy is assumed to be smooth. If 
we allow for perturbations in dark energy then the limits on the 
equation of state parameters changes to —1.25 < w < —0.64 and 
0.20 < Qnr < 0.38 as shown in the top right panel of Figure 5. 

These figures allow us to conclude that: 

• WMAP observations of temperature and polarization 
anisotropies strongly favour models around w = — 1, i.e., the 
ACDM model. As a result, WMAP and Gold+Silver data sets 
have a small region of overlap as the latter does not favour models 
around w = — 1. (It is this disagreement that had led us to suggest 
that the supernova data set could be plagued b y some systematic 
effects dJassal. Bagla. & Padmana bhan 2005a), particularly as it 
contains supernovae from a number of different sources. In that 
work, we have used WMAP first year data.) 

• WMAP and SNLS data sets have a region of overlap within 
68% confidence levels. 
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• There is no significant change in the likelihood contours for 
other cosmological parameters as we go from the cosmological 
constant model to dark energy with a constant equation of state 
parameter (not constrained to w = —1), or when we go from a 
smooth dark energy to the model where dark energy is allowed to 
cluster. 

Thus we can say that SNLS and WMAP data are in (much) better 
agreement as compared to the Gold+Silver and WMAP data sets. 



3.3 Varying w(z) 

It has been claimed that observations, in particular observa- 
tions of high redshift supernovae (the Gold and Gold+Silver 
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20051: IWanell2004 IWang et al.ll2004h. As such a variation is im- 



possible if acceleration of the universe is caused by the cosmo- 
logical constant, it is important to test this claim. Note that the 
term "Evolution of dark energy" has been used for evolution of 
the equation of state parameter, as well as for evolution of energy 
density for the dark energy component. In an earlier study using 
the Gold+Silver data set, we had found that supernova observa- 
tions do not favour evolution of the equation of state parameter over 
models with a constant w -C — 1. But these models are favoured 
strongly as compared to the cosmological constant model, which 
was allowed with V = 6.3% amongst models with constant w. 
When combined with WMAP and other constraints, the allowed 
variation of dark energy is restricted to a narrow range and models 
around the Cosmological constant are favoured (Selia ketalj200l ; 
Jassal, Bagla, & Padmanabhan 2005a). We should note that if we 
combine only the Gold+Silver (or Gold) supernova and WMAP 
data then results favour evolution of pde, but adding observations 
of galaxy clustering removes this inclination. 

We studied constraints on models of varying dark energy with 
the SNLS and WMAP5 data and the results are summarised in ta- 
ble 2, which gives the ranges of parameters allowed at 95% con- 
fidence level. The supernova data constrains wq but does not ef- 
fectively constrain iwi. On the other hand, CMB data constrains an 
effective equation of state and hence indirectly provides constraints 
on w'(z — 0). This is evident from lower panels of Figure 5, where 
we have plotted confidence contours in Q mr — W plane. The CMB 
contours are significantly narrower than those given by the Super- 
nova data. We do not find any significant changes in the likelihood 
contours for other parameters such as n, Q,Nnh 2 , Q.Bh 2 and r. For 
instance the range of QnrIi 2 is given by 0.21 — 0.33 is valid for 
all the models considered here. This is illustrated in Table 2. 

Given that CMB observations constrain only one number, we 
expect that the constraint on models with varying w(z) should 
constrain only w e ff as defined in Eqn.©. We confirm this by 
plotting likelihood contours for models with varying w(z) on the 
Sljvfl — w eff plane. We also plot contours for constant w models 
on the same plane for reference. Figure 6 shows these contours for 
the models with p = 1 and p — 2. We find very strong coinci- 
dence in the contours for models with variable w(z) when plotted 
with w e ff with the contours for models with constant w thereby 
confirming our conjecture. This also provides an easy approach to 
constraining models with variable w(z) without detailed calcula- 



tions, all that one needs is to check whether w e ff is in the range 
allowed by CMB observations for w in the constant w model. 

We return to the issue of variation of dark energy density al- 
lowed by observations. A pictorial representation of results is given 
in Figure 7, where we have plotted Pde(z)/ Pde(z — 0) as a func- 
tion of redshift. Different panels show the evolution of this quantity 
as allowed by the SNLS data set, WMAP5 observations of tem- 
perature and polarization anisotropies in the CMB and combined 
constraints from WMAP5+SNLS. These are plotted for constant w 
(with and without dark energy perturbations), and for variable w 
with p = 1 and p = 2. Dark energy was assumed to be homo- 
geneous in all cases except for the second row that corresponds to 
constant w models with perturbations in dark energy. We can con- 
clude that: 

• Supernova observations are a tight constraint for models with 
constant w, but these are not as strong as CMB constraints. 

• SNLS+WMAP5 data offers tighter constraints than either data 
set and the cosmological constant is allowed with a high probabil- 
ity. This follows from the complementary nature of the two con- 
straints as seen from the orientation of the likelihood contours (e.g., 
see Figure 5). 

• Supernova observations do not constrain evolution of dark en- 
ergy density in models with a variable w. Very large variation in 
dark energy density is allowed by these observations. 

• WMAP5 observations are, in contrast, a much tighter con- 
straint and do not allow significant variation in dark energy. Indeed, 
the variation in dark energy density allowed by WMAP5 observa- 
tions for models with variable w is not significantly larger than that 
allowed for constant w models. 

• We demonstrate that the constraints on dark energy parameters 
for varying w models are the same as the constraints on constant w 
models if we consider w e f / for varying dark energy models. 

• SNLS+WMAP constraints are essentially dominated by the 
WMAP data and follow the same pattern. SNLS observations add 
to the overall constraint by limiting the range of values allowed for 
w . 

We would like to add a note of caution that the analy- 
sis for varying w models does not take perturbations in dark 
energy into account. However, these are more important for 
w S> —1 or models with rapidly va r ying w |Bean&Dorij 
2004; Caldwell, Dave & Steinhardj 1 19981 ; Iweller & Lewis! I200I : 
Hannestad 2005) and such models are not allowed by observational 
constraints. 



4 DISCUSSION 

In this work we studied the SNLS data set and compared the con- 
straints obtained with constraints from WMAP five year data on 
temperature anisotropies in the CMB. We find that the parame- 
ter values favoured by the two data sets have significant over- 
lap and the two sets can be combined to put tight constraints 
on models of dark energy. In an earlier work we had noted that 
the Gold+Silver data set does not agree with WMAP observa- 
tions in that these favour distinct pa rts of the parameter space 
dJassal, Bagla. & Padmanabhanl l2005a). Constraints from WMAP 
and structure formation favour similar models, but ones distinct 
from those favoured by Gold+Silver supernova observations. This 
indicates some degree of inconsistency between the supernova and 
other observations and it led us to suggest that the Gold+Silver 
data set may be affected by as yet unknown systematic errors 
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Figure 6. The black contours in this figure show constraints on the effective dark energy equation of state w e e / for p = 1 and p = 2 from the WMAP5 data. 
The green lines are constraints on constant w — Cjvjj (without perturbations). These clearly show that the CMB data constrains effective equation of state at 
the last scattering surface. 



Table 2. This table lists the range of parameters allowed within 95% confidence limit from SNLS, WMAP1, WMAP3 and SNLS+WMAP3. 



Parameter 


ACDM 


w=const. 


w=const. 


P= 


1 


P= 


2 














with perturbations 


















-1.92- 


--0.74 


-1.92- 


—0.74 


-1.89- 


—0.61 


-1.9— 


-0.59 


SNLS 


w 






-1.39- 


- -0.58 


-1.63- 


--0.66 


-1.64- 


--0.42 


-1.93- 


- -0.43 


WMAP1 








-1.25- 


- -0.7 


-1.25- 


--0.64 


-1.62 — 


- -0.44 


-1.62 - 


- -0.43 


WMAP5 








-1.47- 


- -0.83 


-1.57- 


--0.88 


-1.46- 


--0.81 


-1.74- 


--0.77 


SNLS+WMAP1 








-1.1- 


—0.9 


-1.1- 


--0.9 


-1.3 — 


--0.8 


-1.42- 


--0.66 


SNLS+WMAP5 
















-4.82 


—3.3 


-4.79- 


—4.23 


SNLS 


w'(z = 0) 














-3.09- 


-1.32 


-2.5- 


-4.87 


WMAP1 
















-1.8- 


-1.2 


-3.3- 


-2.7 


WMAP5 
















-0.99- 


-1.04 


-2.22- 


-4.79 


SNLS+WMAP1 
















-1.25- 


-0.91 


-2.6- 


-2.7 


SNLS+WMAP5 




0.22- 


-0.31 


0.11- 


-0.47 


0.11- 


-0.47 


0.11- 


-0.48 


0.11- 


-0.48 


SNLS 




0.20- 


-0.45 


0.16- 


-0.43 


0.20- 


-0.47 


0.17- 


-0.45 


0.18- 


-0.44 


WMAP1 




0.19- 


-0.33 


0.19- 


-0.38 


0.19- 


-0.39 


0.19- 


-0.39 


0.19- 


-0.38 


WMAP3 




0.21- 


-0.34 


0.2 — 


0.38 


0.2 — 


0.38 


0.2 — 


0.38 


0.20- 


-0.38 


WMAP5 




0.22- 


-0.31 


0.15- 


-0.36 


0.18- 


-0.41 


0.16- 


-0.38 


0.19- 


-0.39 


SNLS+WMAP1 




0.22 


— 0.3 


0.22- 


-0.3 


0.22- 


-0.3 


0.22- 


-0.33 


0.21- 


-0.32 


SNLS+WMAP5 


h 


0.61- 


-0.79 


0.61- 


-0.78 


0.6 — 


0.79 


0.6 — 


0.78 


0.61- 


-0.78 


WMAP1 




0.66- 


-0.77 


0.60- 


-0.78 


0.60- 


-0.78 


0.61- 


-0.79 


0.61- 


-0.78 


WMAP5 




0.69- 


-0.77 


0.68- 


-0.78 


0.68- 


-0.79 


0.67- 


-0.77 


0.65- 


-0.78 


SNLS+WMAP1 




0.68- 


-0.74 


0.68- 


-0.74 


0.68- 


-0.74 


0.66- 


-0.76 


0.65- 


-0.76 


SNLS+WMAP5 


n B h 2 


0.02- 


-0.027 


0.021- 


-0.028 


0.02 — 


0.027 


0.02 — 


0.027 


0.02 — 


0.027 


WMAP1 




0.021- 


-0.024 


0.021- 


-0.023 


0.02 — 


0.024 


0.021 — 


0.0235 


0.022- 


-0.024 


WMAP5 




0.02- 


-0.027 


0.021- 


-0.027 


0.021- 


-0.027 


0.021- 


-0.027 


0.021- 


-0.028 


SNLS+WMAP1 




0.021- 


-0.023 


0.021- 


-0.024 


0.021- 


-0.024 


0.021- 


-0.023 


0.021- 


-0.024 


SNLS+WMAP5 


n. 


0.93- 


-1.08 


0.93- 


-1.1 


0.93- 


-1.09 


0.93- 


- 1.1 


0.93 — 


1.098 


WMAP1 




0.94- 


-0.99 


0.93- 


-0.99 


0.93- 


-0.99 


0.93- 


-0.99 


0.93- 


-0.99 


WMAP5 




0.93- 


-1.09 


0.94- 


-1.08 


0.94- 


-1.09 


0.93- 


-1.09 


0.93 — 


1.097 


SNLS+WMAP1 




0.94- 


-0.99 


0.94- 


-0.99 


0.93- 


-0.99 


0.93- 


-0.99 


0.94- 


-0.99 


SNLS+WMAP5 


T 


0.002 


— 0.33 


0.011- 


-0.39 


0.007- 


-0.35 


0.13- 


-0.4 


0.016- 


-0.39 


WMAP1 




0.054 


— 0.12 


0.05- 


-0.12 


0.055 - 


-0.12 


0.054- 


-0.12 


0.054 - 


-0.12 


WMAP5 




0.004 


— 0.33 


0.008- 


-0.34 


0.045- 


-0.37 


0.013- 


-0.38 


0.017- 


-0.396 


SNLS+WMAP1 




0.054 


— 0.12 


0.053- 


-0.12 


0.055- 


-0.13 


0.52- 


-0.12 


0.05 — 


0.125 


SNLS+WMAP5 
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Figure 7. This figure shows allowed range in variation of dark energy density as a function of redshift. The top row is for homogeneous dark energy model, 
the second row is for perturbed dark energy, the third row is for varying w with p = 1, and, the last row is for varying w with p = 2. The white region is 
the allowed range at 63% confidence level, the hatched region is the one disallowed range of dark energy density at 95% confidence level and the solid (blue) 
region is the one ruled out at 99% confidence level. In all the rows, the left most plot shows range allowed by SNLS data, the middle one with WMAP5 data 
and the right one shows the range allowed by combined data. 
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jjassal. Bagla. & Padmanabhanl 120053) . One reason for doubting 
the supernova data is the heterogeneit y of sources from which 
the particular dat a set was collected dRiess et al . 2004). SNLS 
dAstier et alj2005l) is a homogeneous data set and should not suffer 
from such problems and indeed we find that there is no inconsis- 
tency between SNLS and WMAP observations. 

This highlights the usefulness of CMB observations for 
constraining models of dark energ y jEisenstein & White! [20041 ; 
Jassal, Bag la & Padma nabhan 2005J). We believe that CMB obser- 
vations should be used for testing any model of dark energy as su- 
pernova observations do not constrain models with varying w ef- 
fectively0. Thus one should use CMB observations as well and not 
rely only on supernova observations for constraining such models. 

We would like to add a note of caution against combining the 
SNLS data with other data sets of high redshift supernovae in light 
of the very different nature of these data sets. Indeed, one should 
use homogeneous data sets like the SNLS in isolation to avoid the 
problems mentioned above. 

In terms of models, we find that the cosmological constant 
is favoured by individual observations (SNLS and WMAP) as 
well as in the combined data set with very high probability. 
Table 2 gives allowed values of all cosmological parameters at 
95% confidence level by SNLS, WMAP as well as the combined 
data set. For the cases where a similar analysis has b een done by 
other s , our results are con s istent with other findings i Bridle et al. 
20031 ; ISpergel et all 120031; iMaccio et alj [20031: iLinder & Jenkins! 
2003; IPogosvan Bond & Contaldil 120031; iTegmark et all 1 2004 ; 
Wang & Tegmarkl |2004 iGiovi, Baccigalupi & Perrottal 12004 ; 
Hannestad 1 |2004 iHuterer & Cooravl 12004 iLed 120051; iLee & 



2003 k ILee. Lee & Ng 
20051; iRapetti. Allen & Welle. 



20031; 



Nesseris & PerivolaropoulosI 2005|; Seliaket al] 
Shen. Wang. Abdalla & Sul 120051; iFeng, Wang & Zhand 
Xia et alj I2OO6 Jbt lAmendola. Campos. & Rosenfeld 



lettol 



Mainini, Colombo & Bonomett* 

|2004 IPogosvarj ~ 
20051; ISeliaketalJ 



Calvo & Maroto 12006); ICa rneiro et all 120061; iDantas et al] 



2004 



2005 



2005 



2007 



2007 



Elizalde et al.l 1200 8l : iMota. Kristiansen. Koivisto & Groenebooml 
20071) . 

We have discussed the origin of the constraint on dark energy 
models from CMB observations at length. We may conclude from 
the analysis presented here that: 

• Location of acoustic peaks in the angular power spectrum of 
CMB anisotropies is the main source of constraints. 

• CMB observations only constrain w e ff, an effective value of 
the equation of state parameter defined in Eqn.{6}. This can be used 
to translate constraints on models with constant w to models where 
dark energy properties vary with time. 

• We have discussed models with fin = 1 in this paper. In 
case this constraint is relaxed then the well known degeneracy 
between w and fin loosens the constraints. However, it is well 
known that the SN and CMB data are complementary and can 
be combined to provide fairly tight constraints even in t his case 



dPerlmutter. Turner & White! 1999f ; lHuterer & Turnerl200lh . We do 
not expect variations in curvature to modify our conclusions about 
W e ff being the only dark energy related quantity constrained by 
CMB observations. 

• Integrate Sachs-Wolfe (ISW) effect due to perturbations in 
dark energy can, in principle, lead to variations in the CMBR angu- 
lar power slectrum at small I. Our analysis of models with variable 
w without perturbations does not take this into account. Various 

3 The main constraint on varying w models is from the CMB data 



analyses have shown that ISW does not contribute significantly to 
constrainin g cosmologi cal parameters from CMB data (see for ex- 
ample, dXia et ai] |2009h ). The main reason for this is that ISW af- 
fects power at small I, and due to large cosmic variance these modes 
do not contribute much to the overall likelihood. 



ACKNOWLEDGEMENTS 

Numerical work for this study was carried out at cluster 
computing facility in the Harish-Chandra Research Institute 
(http://cluster.hri.res.inl. This research has made use of NASA's 
Astrophysics Data System. 



REFERENCES 

Andrianov A. A., Cannata F, Kamenshchik A. Y., 2005, Phys. 

Rev. D 72, 043531 
Aguirregabiria J. M., Lazkoz R., hep-th/0402190 
Alam U„ Sahni V., Starobinsky A. A., 2004, JCAP, 0406, 008 
Alam U., Saini T.D., Sahni V., Starobinsky A. A., 2004, Mon. Not. 

Roy. Astron. Soc, 354, 275 
Alcaniz, J. S., 2006, Braz. J. Phys., 36, 1109 
Amendola L., Quercellini C, 2003, Phys. Rev. D 68, 023514 
Amendola L., Campos G. C, Rosenfeld R., 2007, PhRvD, 75, 

083506 

Apostolopoulos P. S., Tetradis N., 2006, PhRvD, 74, 064021 
Armendariz-Picon C, Mukhanov V., Steinhardt P. J., 2001, 

PhRvD, 63, 103510 
Arianto, Zen F. P., Gunara B. E., Triyanta, Supardi, 2007, JHEP, 

9,48 

Astier P. et al. 2005, |astro-ph/05 1 0447 1 

Bagla J. S., Padmanabhan T, Narlikar J. V., 1996, Comments on 

Astrophysics 18, 275; |astro-ph/9511102| 
Bartolo N., Corasaniti P.-S., Liddle A. R., Malquarti M., 2004, 

PhRvD, 70, 043532 
Jassal H. K., in preparation 

Bagla J. S., Jassal H. K., Padmanabhan T, 2003, Phys. Rev. D 

67, 063504; |astro-ph/0212198| 
Barris B. J., et al., 2004, ApJ, 602, 571 
Bassett B. A., |astro-ph/0407201| 

Bassett B. A., Corasaniti P. S., Kunz M., 2004, Astrophys. J., 617, 
LI 

Bean R., Dore O., 2004, Phys. Rev. D 69, 083503 
Benabed K., Bernardeau F, 2001, Phys. Rev. D 64, 083501 
Bento M. C, Bertolami O., Sen A. A., 2002, Phys. Rev. D 66, 
043507 

Bludman S. A., Roos M., 2002, Phys. Rev. D 65, 043503 
Bridle S. L., Lahav O., Ostriker J. P., Steinhardt P. J., 2003, Sci, 
299, 1532 

Briscese F, Elizalde E., Nojiri S., Odintsov S. D., 2007, PhLB, 
646, 105 

Bronnikov K. A., Fabris J. C, Goncalves S. V. B., 2007, JPhA, 
40, 6835 

Bronnikov K. A., Starobinsky A. A., 2007, JETPL, 85, 1 
Burgess C. P., 2003, Int. J. Mod. Phys. D 12, 1737 
Calvo G. B., Maroto A. L., 2006, PhRvD, 74, 083519 
Carneiro S., Pigozzo C, Borges H. A., Alcaniz J. S., 2006, 

PhRvD, 74, 023532 
Carroll S. M., Press W. H., Turner E. L., 1992, ARA&A, 30, 499 



CMB constraints on Dark Energy 1 1 



Carroll S. M, Hoffman M, Trodden M., 2003, Phys. Rev. D 68, 
023509 

Caldwell R. R., 2002, Phys. Lett. B 545, 23 
Caldwell R. R., Dave R., Steinhardt P. J., 1998, Phys. Rev. Lett. 
80, 1582 

Chiba T., 2002, Phys. Rev. D 66, 063514 

Chimento L. P., Feinstein A., 2004, Mod. Phys. Lett. A 19, 761 

Choudhury T. R., Padmanabhan T.,Astron.Astrophys., 429: 807, 

(2005) |astro-ph/03ll622) 
Choudhury D , Ghoshal D., Jatkar D. P., Panda S., 2002, Phys. 

Lett. B 544, 231 
Clarkson C, Cortes M., Bassett B., 2007, JCAP, 8, 11 
Cline J. M., Jeon S. Y, Moore G. D., 2004, Phys. Rev. D 70, 

043543 

Cognola G., Elizalde E., Nojiri S., Odinstov S. D., Zerbini S., 
|hep-th/0601008l 

Copeland, E. J. and Sami, M. and Tsujikawa, S., 2006, Int. J. Mod. 

Phys. D 15, 1753 
Corasaniti P. S., Kunz M., Parkinson D., Copeland E. J., Bassett 

B. A., |astro-ph/0406018| 
Dabrowski M. P., Stachowiak T., Szydlowski M., 2003, Phys. Rev. 

D68, 103519 
Daly R. A., Djorgovski S. G., |astro-ph/0405063| 
Dantas M. A., Alcaniz J. S., Jain D., Dev A., 2007, A&A, 467, 

421 

Dunkley J. et al., Astrophys. J. Suppl. 180, 306 

Das S.,Banerjee N, Dadhich N., |astro-p h/0505096 

Dev A., Jain D., Alcaniz J. S., 2003, Phys. Rv. D 67, 023515 

Doran M., Robbers G., 2006, JCAP, 6, 26 

Efstathiou G., Sutherland W. J., Maddox S. J., 1990, Nature 348, 
705 

Eisenstein D. J., White M., 2004, PhRvD, 70, 103523 

Elizalde E., Nojiri S., Odintsov S. D„ 2004, PhRvD, 70, 043539 

Elizalde E., Nojiri S., Odintsov S. D., Saez-Gomez D., Faraoni V., 

2008, arXiv, 803, larXiv:0803.1311| 
Ellis J. R., 2003, Phil. Trans. Roy. Soc. Lond. A 361, 2607 
Fabris J. C, Shapiro I. L., Sola J., 2007, JCAP, 2, 16 
Fabris J. C, Goncalves S. V. B., 2006, PhRvD, 74, 027301 
Feng B., Wang X., Zhang X., 2005, Phys. Lett. B 607, 35 
Frampton P. H., |hep-tn/0302007| 

Frolov A. V., Kofman L., Starobinsky A. A., 2002, Phys. Lett. B 
545, 8 

Garnavich P. M., et al., 1998, ApJ, 509, 74 

Gibbons G. W., 2002, Phys. Lett. B 537, 1 

Gibbons G. W., 2003, Class. Quant. Grav. 20, S321 

Gibbons G. W., |hep-th/0302199| 

Giovi F, Baccigalupi C, Perrotta F. , [astro-ph/04 1 1 702 1 

Gong Y. G., |astro-ph/0401207| 

Gonzalez-Diaz P. F, 2002, Phys. Rev. D 62, 023513 

Gonzalez-Diaz P. F, 2000, Phys. Lett. B 481, 353 

Gonzalez-Diaz P. F, 2003, Phys. Rev. D 68, 021303 

Gorini V., Kamenshchik A. Y, Moschella U., Pasquier V., 2004, 

PhRvD, 69, 123512 
Gordon C, Hu W„ 2004, PhRvD, 70, 083003 
Gordon C, Wands D., 2005, PhRvD, 71, 123505 
Gorini V., Kamenshchik A. Y, Moschella U., Pasquier V., 2001, 

Phys. lett. B 511, 265 
Hannestad S., 2004, JCAP, 0409, 001 
Hannestad S., 2005, Phys. Rev. D 71, 103519 
Hannestad S., Mortsell E., 2004 JCAP, 0409, 001 
Hao J. g. Li X. z., 2003, Phys. Rev. D 67, 107303 
Holman R. and Naidu S. 2004, astro-ph/0408102 



Hu W., 2005, PhRvD, 71, 047301 

Hu W., Dodelson S., 2002, ARA&A, 40, 171 

Huterer D., Turner M. S., 2001, PhRvD, 64, 123527 

Huterer D., Cooray A., astro-ph/0404062 

Jain P. and Ralston J. P. 2005, ApJ In Press, |astro-ph/0506478| 

Jassal H. K., 2004, Pramana 62, 757 

Jassal H. K., |hep-th/03 122531 

Jassal H. K., Bagla J. S., Padmanabhan T, 2005, Mon. Not. Roy. 

Astron. Soc. 356, Lll-L16; |astro-ph/ 0404378 
Jassal H. K., Bagla J. S., Padmanabhan T, 2005, PhRvD, 72, 

103503; astro-ph/0506748 
Jassal H. K., |arXiv:0903.5370l 
Jimenez R., 2003, New Astron. Rev. 47, 761 
Jimenez R., Verde L., Treu T, Stern D., 2003, Astrophys. J. 593, 

622 

Jonsson J., Goobar A., Amanullah R., Bergstrom L., 
|astro-ph/0404468| 

Komatsu E., et al., 2009, ApJS, 180, 330 

Kowalski M., et al., 2008, ApJ, 686, 749 

Kim C. j., Kim H. B., Kim Y. b., 2003, Phys. Lett. B 552, 111 

Lazkoz R., Nesseris S., Perivolaropoulos L., astro-ph/0503230 

Lee S., |astro-ph/0504650| 

Lee W, Ng K., 2003, Phys. Rev. D 67 107302 

Lee D., Lee W., Ng K., |astro-ph/0309316| 

Li M., 2004, Phys. Lett. B603 

Li H., Xia J., Zhao G, Fan Z., Zhang X., 2008, ApJL, LI 
Linder E. V., Jenkins A., 2003, Mon. Not. Roy. Astron. Soc. 346, 
573 

Linder E. V., Robbers G., 2008, JCAP, 0806, 004 
de la Macorra A., Piccinelli G., 2000,Phys. Rev. D 61, 123503 
Maccio A. V., Bonometto S. A., Mainini R., Klypin A., 
|astro-ph/0309439| 

Mainini R., Colombo PL., Bonometto S. A., astro-ph/05030361 
Mainini R., 2009, JCAP 0904, 017 
Mainini R., 2008, JCAP 0807, 003 

Malquarti M., Copeland E. J., Liddle A. R., Trodden M., 2003, 

Phys. Rev. D 67, 123503 
Melchiorri A. et al. [Boomerang Collaboration], 2000, Astrophys. 

J. 536, L63 
Milton K. A., 2003, Grav. Cosmol. 9, 66 
Mota D. F, van de Bruck C, 2004, A&A, 421, 71 
Mota D. F, Kristiansen J. R., Koivisto T, Groeneboom N. E., 

2007, Mon. Not. R. Astron. Soc. 382, 793 
Nesseris S., Perivolaropoulos L., 2005, Phys. Rev. D 72, 123519 
Nojiri S„ Odintsov S. D., 2003, Phys. Lett. B 562, 147 
Nojiri S., Odintsov S. D., Tsujikawa S., 2005, PhRvD, 71, 063004 
Nunes N. J., Mota D. F, 2006, MNRAS, 368, 751 
Onemli V. K ., Woodard R. P., 2002, Class. Quant. Grav. 19, 4607 
Onemli V. K., Woodard R. P., 2004, Phys. Rev. D 70, 107301 
Ostriker J. P., Steinhardt P. J., 1995, Nature 377, 600 
Padmanabhan T, 2002, PhRvD, 66, 021301; [hep^th/0204150| 
Padmanabhan T, 2002, Class. Quant. Grav. 19, 5387; 

|gr-qc/02040T9l 

Padmanabhan T, 2003, Phys. Rept. 380, 235; |hep-th/0212290| 
Padmanabhan T, Choudhury T. R., 2003, Mon. Not. Roy. Astron. 

Soc. 344, 823; |astro-ph/0212573| 
Padmanabhan T, 2005, Curr. Sci. 88, 1057; astro-ph/04 110441 
Padmanabhan T, Understanding Our Universe: Current Status 

and Open Issues; gr-qc/0503107 
Padmanabhan T, 2005, Phys. Rept. 49, 406; |gr-qc/0311036| 
Padmanabhan T, astro-ph/05 10492 
Padmanabhan T, gr-qc/0510015 



12 Jassal, Bagla and Padmanabhan 



Peebles P. J. E., Ratra B., 2003, Rev. Mod. Phys. 75, 559 
Peiris H. V. and Spergel D. N. 2000, ApJ 540, 605 
Perivolaropoulos L., astro-ph/0601014 
Percival W. J., et al., 2007, ApJ, 657, 645 
Perlmutter S., et al., 1999, ApJ, 517, 565 

Perlmutter S., Turner M. S., White M., 1999, Phys. Rev. Lett. 83, 

670 

Perlmutter S. and Schmidt B. P. 2003, astro-ph/0303428 
Pogosyan L., astro-ph/0409 059| 

Pogosyan D., Bond J. R., Contaldi C. R., |astro-ph/0301310| 

Polarski D., Ranquet A., 2005, Phys. Lett. B 627, 1 

Rapetti D., Allen S. W„ Weller J.,jastro-ph/0409574l 

Ren J., Meng X., |astro-ph/0511163| 

de Ritis R., Marino A. A., 2001, Phys. Rev. D 64, 083509 

Riess A. G., et al., 1998, AJ, 116, 1009 

Riess A. G., et al., 2004, ApJ, 607, 665 

Rubano C, Scudellaro P., 2002, Gen. Rel. Grav. 34, 307 

Sahni V., Starobinsky A. A., 2000, Int. J. Mod. Phys. D 9, 373 

Sahni V., Shtanov Y, 2003, JCAP 0311, 014 

Sahni V., Shafieloo A., Starobinsky A., 2008, Phys. Rev. D78, 

103502 

Scherrer R. J., |astro-ph/0402316| 
Seljak U. and Zaldarriaga M. 1996, ApJ 469, 437. 
Seljak U., et al., 2005, PhRvD, 71, 103515 
Sen A., |hep-th/0312T53l 

Sen S., Seshadri T. R., 2003, Int. J. Mod. Phys. D 12, 445 
Sen A. A., Scherrer R. J., 2005, Phys. Rev. D 72, 0635 1 1 
Sergijenko O., Kulinich Y., Novosyadlyj B., Pelykh V., 2008, 

arXiv, |arXiv:0809.3349l 
Gupta S., Saini T. D., Laskar T., 2007, [astro-ph/0701683| 
Shen J., Wang B., Abdalla E., Su R., 2005, Phys. Lett. B 609, 200 
Shiu G., Wasserman I., 2002, Phys. Lett. B 541, 6 
Singh P., Sami M., Dadhich N., 2003, Phys. Rev. D 68, 023522 
Sola J., Stefancic H., |astro-ph/0507110| 
Spergel D. N. et al., 2003, ApJS 148, 175 

Steinhardt P. J., Royal Society of London Philosophical Transac- 
tions Series A, 361, 2497 
Subramanian K., |astro-ph/0411049| 
Tegmark M., et al., 2004, ApJ, 606, 702 
Tonry J. L. et al., 2003, ApJ 594, 1 

Unnikrishnan S., Jassal H. K., Seshadri T. R., 2008, 
larXiv:0801. 20171 

Urena-Lopez L. A., Matos T, 2000, Phys. Rev. D 62, 081302 

Uzawa K., Soda J., 2001, Mod. Phys. Lett. A 16, 1089 

Virey J.-M., Taxil P., Tilquin A., Ealet A., Tao C, Fouchez D., 

2005, PhRvD, 72, 061302 
Wang Y., |astro-ph/0404484| 

Wang Y, Tegmark M., 2004, Phys.Rev.Lett.92, 241302 
Wang Y. et al., |astro-ph/0402080| 
Weinberg S. 1989, Rev.Mod.Phys. 61, 1 

Weller J., Lewis A. M., Mon. Not. R. Astron. Soc, 2003, 346, 

987 

White M., Cohn J. D., [astro-ph/0203120| 

Wood-Vasey W. M., et al., 2007, ApJ, 666, 694 

Xia J., Zhao G., Feng B., Li H., Zhang X., 2006, Phys. Rev. D 73, 

063521 

Xia J., Zhao G., Li H., Feng B., Zhang X., astro-ph/0605366 
Xia J., Li H„ Zhao G., Zhang X., 2008, Phys. Rev. D78, 0835254 
Xia J.-Q., Viel M., Baccigalupi C, Matarrese S., 2009, JCAP, 9, 3 



